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ABSTRACT 

Context. Molecular hydrogen observations towards Herbig-Haro objects provide the possibility of studying physical processes related to star 
formation. 

Aims. Observations towards the luminous IRAS source IRAS 11101-5928 and the associated Herbig-Haro objects HH135/HH136 are obtained 

to understand whether high-mass stars form via the same physical processes as their low-mass counterparts. 
' Methods. Near-infrared imaging and spectroscopy are used to infer H2 excitation characteristics. A theoretical H2 spectrum is constructed from 

a thermal ro- vibrational population distribution and compared to the observations. 
, Results. The observations reveal the presence of a well-collimated, parsec-sized Ht outflow with a total H2 luminosity of about 2Lo. The bulk 

of the molecular gas is characterized by a ro-vibrational excitation temperature of 2000 + 200 K. A small fraction (0.3%) of the molecular gas 
] is very hot, with excitation temperatures around 5500 K. The molecular emission is associated with strong [Fell] emission. The Ht and [Fell] 
. emission characteristics indicate the presence of fast, dissociative J-shocks at speeds of Vs ~ 100 km s"'. Electron densities of = 3500-4000 

cm"' are inferred from the [Fell] line ratios. 

Conclusions. The large H2 luminosity combined with the very large source luminosity suggests that the high-mass protostar that powers the 
HH135/HH136 flow forms via accretion, but with a significantly increased accretion rate compared to that of low-mass protostars. 

Key words. ISM:individual objects: HH135/HH136 - ISM: Herbig-Haro objects - ISM: jets and outflows 



1. Introduction 

Near-infrared studies of molecular hydrogen emission in 
Herbig-Haro (HH) objects provide a powerful tool to gain 
insight in the physical processes that occur during the early 
phases of low-mass star formation. The total H2 luminosities 
are proportional to the accretion rates in the early phases of the 
protostellar evolution (Stanke l20001 Froebrich et al. 2003 1, and 
the H2 ro-vibrational population distribution allows us to deter- 
mine the physical scenarios that are at work in the supersonic 
jets that form during the accretion phase (McCoev 120041 and 
references therein). A complementary tool to study such re- 
gions is available via the near-infrared emission lines of [Fell] 
and other atomic emission lines, that allow us to measure elec- 
tron densities in the shocked material (e.g., Nisini et al. 2002). 

In a number of HH objects, and most prominently in HH 
objects that show pronounced [Fell] emission, a temperature 
stratification is inferred from the H2 observations, where part of 
the molecular gas reach temperatures above 5000 K (Giannini 

Send offprint requests to: R. Gredel 

* Based on observations collected at the European Southern 
Observatory, La Silla, Chile. Figures .1-.4 and Tables .1-.9 are only 
available in electronic form via _http://www.edpsciences.org_ 



et al. 12004 1. The combined H2 and [Fell] emission has been 
explained in terms of fast J-type shocks, where the [Fell] emis- 
sion is produced in dissociative parts of the shocks, and where 
H2 arises in the slower, non-dissociative regions (e.g., Gredel 
1994i. Detailed model calculations of J-type shocks with mag- 
netic precursors confirm such a model, and produce popula- 
tion distributions where the rotational excitation temperatures 
increase with increasing vibrational state of H2 (e.g.. Flower 
et al. 2003 1. C-type shocks at largely different physical con- 
ditions produce similar H2 population distributions, however, 
and the H2 population distribution alone does not allow us to 
distinguish between both scenarios (Flower et al. I2003I I. The 
presence of [Fell] emission is generally interpreted in terms of 
dissociative J-type shocks, although it is conceivable that C- 
type shocks produce [Fell] as well (Le Bourlot et al. l2002> . 

Near-infrared studies of outflows from intermediate- and 
high-mass star forming regions are rare, either because such 
outflows are not very frequent per se or because high-mass 
star-forming regions are deeply embedded in general. Some 
outflows that are observed from high-luminous sources such 
as Orion (IO^'Lq) lack the high degree of collimation that 
is typical for outflows from low-mass star-forming regions. 
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Other flows that emerge from luminous IRAS sources, such as 
HH80/81 (21O'*L0, Marti, Rodriguez, & Reipurth [l995t . IRAS 
1 6547-4247 (6 1 (fl^. Brooks et al. 2003 ), or IRAS 1 8 1 5 1 - 1 208 
(210"* Lq, Davis et al. 2004), are highly collimated. It is not clear 
whether intermediate- and high-mass stars form in a similar 
way to low-mass stars, but with enhanced accretion rates, or 
whether different processes, such as the merging of low-mass 
protostars, are at work. Enlarging the sample of near-infrared 
observations of regions of intermediate- and high-mass star for- 
mation and comparing the general properties of their H2 and 
[Fell] emission with those of low-mass star-forming regions is 
therefore desirable. 

In the following, a study of the molecular outflow 
HH135/HH136, which is powered by the cold and very lu- 
minous (lO'^Lg) IRAS source IRASl 1 101-5928, is presented. 
The pair of Herbig-Haro objects HH135/HH136 was discov- 
ered by Ogura & Walsh S\992\ in an objective prism survey 
and is located in the bright rimmed cloud No. 64 of Sugitani 
& Ogura J1994> in the Eastern Carina region. The HII re- 
gion, also known as Gum 36, is believed to be excited by the 
open cluster Stock 13, for which photometric distances of 2.7 
kpc are available (Steppe 1977). The general morphology of 
HH135/HH136 indicates that the two objects are formed at the 
opposite flow directions of a bipolar flow, which is driven by 
IRASl 1 101-5829. The velocity field of the HH135/HH136 re- 
gion is complex, however, and the fact that the main part of 
the emission from both HH135 and HH136 is blue shifted led 
Ogura & Walsh (1992) to conclude that HH135 and HH136 
form two different, independent flows. Near-infrared JHK po- 
larimetric observations of the associated reflection nebula car- 
ried out by Tamura et al. ( flWTb showed that IRAS 1 1 101-5829 
is the only illuminating source of the nebula that is associ- 
ated with HH135/HH136. A more recent millimetre study by 
Ogura et al. (1998 1 presented a model that explains the ob- 
served velocity features and where HH135 and HH136 form 
part of a single, bipolar flow driven by IRASl 1 101-5829. This 
view is supported by the very recent polarimetric observations 
by Chrysostomou et al. (2005 1, who proposed that a strong he- 
lical magnetic field threading though HH135/HH136 maintains 
the strong collimation of the flow. 

The purpose of this paper is to investigate the molecular 
hydrogen emission of HH135/HH136, and to study how the 
emission characteristics differ from those of low-mass star- 
forming regions. The observations are presented in Sect. |2l 
which also includes a more detailed discussion of the data re- 
duction method and errors. The imaging results are given in 
Sect. 13.11 and the spectroscopy in Sect. 13.21 Section 13.31 con- 
tains an analysis of the global properties of the molecular hy- 
drogen emission, and Sect. 13.41 summarises the results of an 
analysis derived from the [Fell] lines, which are observed to- 
wards HH135/HH136. The main conclusions of this study are 
summarized in Sect.|3] 

2. Observations and reduction 

Near-infrared imaging and spectroscopy of the HH135/HH136 
complex was carried out in March 1999 and February 2004 at 
the La Silla Observatory, using SOFI at the New Technology 
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Fig.l. [SII] image of HH135/HH136 and part of the G 36 
region in Carina. Optical knots in HH136 are labeled fol- 
lowing the notation by Ogura & Walsh \\992\ . The crossed 
square south of knot J marks the position of the IRAS source 
IRASl 1 101-5829. North is up and east is left. 
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Fig. 2. Gunn z image of HH135/HH136. 

Telescope (NTT). Images were obtained at both epochs in three 
narrow-band filters NB164, NB212, and NB228, using a pixel 
scale of O'.'144/px. Central wavelengths are 1.644 jum, 2.124 
jum, and 2.280 jim, respectively, and widths are 0.025 /zm, 
0.028 jum, and 0.030 yum, respectively, for the NB164, NB2I2, 
and NB228 filters, according to the information given in the 
ESO web pages. Dithered images were used to construct back- 
ground subtracted images, including corrections for sky, dark- 
current, and bias. This method assumes that the background is 
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constant and does not vary across an image, which is not the 
case in the H-band. All images were re-centered to a common 
reference, but no effort was made to correct for field distortion 
across the image. 

In March 1999, spectroscopic observations were obtained 
along three slit positions using the red grism GR in order 1, 
covering the 1.52-2.52 // wavelength region. The slit was set 
to a width of 0'.'6, and measurements of the width of individual 
spectral lines showed that an effective spectral resolution of R 
= AjhA- 550 was achieved. Spectroscopic observations were 
also carried out in February 2004, using the blue grism GB in 
order 1 and a slit of 0'.'6, covering the 0.94-1 .60 yum wavelength 
region at a spectral resolution of R=850. Additional spectra 
were gathered using the HR grism in orders 3 and 4, covering 
the 2.0-2.3 jim and 1.5-1.8 fim wavelength regions, respec- 
tively. The K- and H-band observations obtained with the HR 
grism were carried out with a slit width of 0'.'6 as well. Spectral 
resolutions, as obtained from the width of the emission lines, 
were R=1066 and R=1600, respectively. Explicit sky frames 
were obtained by moving the object along the slit. This is pos- 
sible because the infrared emission of HH135/HH136 is ex- 
tended over about 1 arcmin, compared to the total slit length of 
5 arcmin. 

The spectroscopic observations were reduced using the 
MIDAS long package. The long-slit spectra suffer from a pro- 
nounced curvature of the spectral lines perpendicular to the 
dispersion direction, and the distorted 2-D frames were recti- 
fied by use of a second-order polynomial. Pixel units along the 
dispersion direction were then transformed into a linear wave- 
length scale by the use of Xenon- Argon calibration frames. 
A first-order background subtraction was obtained by build- 
ing (object - sky) difference frames. Short-term temporal 
and spatial variations at scales comparable to the length of 
the slit result in residuals in the background in a large num- 
ber of the difference frames. A second-order sky subtraction 
was achieved via the standard technique of extracting the posi- 
tive and negative spectra from the difference frames and adding 
them in the correct sense. Extraction windows to obtain the 1-D 
spectra from the 2-D difference frames were typically 12 pix- 
els along the slit. On the spatial scale of O'.'29/px, the extraction 
windows correspond to solid angles Q of about 2 arcsec^ on the 
sky or Q = 510"" sr The atmospheric transmission was deter- 
mined via observations of standard stars at slit widths of Qf!6. 
Stellar fluxes were measured with a second integration, but slits 
opened to 2" to minimise slit losses. During the observations, 
the seeing varied from 0'.'5 to r.'5. In phases of poor seeing, 
the flux conversion factors determined from the stellar obser- 
vations varied by up to 30%. The highly peaked nature of the 
line emission in some of the knots, combined with positional 
errors in the slit position for the various grism settings and the 
narrow slit of 0'.'6 adopted here, causes differences in the re- 
ported fluxes towards some of the knots that reach, in the worst 
cases, a factor of about three (cf. Tables Q]and0. This is well 
illustrated in Fig.|8l which shows that an error of less than 0'.'5 
in the position of the slit would cause the loss of most of the 
flux in knot B 1 . The consequences of these kinds of flux un- 
certainties for the analysis of the H2 and atomic emission lines 
in HH135/HH136 are discussed in detail in Sects. I32landl3l1 
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Fig. 3. Narrow-band image of the HH135/HH136 region ob- 
tained in the NB212 filter, which is centered on the (1,0) S(l) 
H2 emission line near 2.12 jim. The three double lines indi- 
cate slit locations and slit widths for the February 2004 spec- 
troscopic observations, with position angles of PA = 33°, 43°, 
and 45°. North is up and east is left. 

The tables of Pickles ("1998^ were used to convert the counts of 
the standard stars to flux units. 

In addition to the near-infrared observations, optical im- 
ages in the [SII] (/Icenter - 673nm, FWHM = 6nm) and 
Gunn z(/lcut-on - 840 nm) were obtained in April 1993, using 
EFOSC2 at the 2.2m telescope and CCD#19. Integration times 
were 300 sec in [SII] and 2100 sec in Gunn z. 

3. Results and discussion 

The narrow-band [SII] image of the HH135/HH136 region and 
the Gunn z image are shown in Figs.^andl^ The observations 
were obtained during an atmospheric seeing of about r.'3. The 
image was re-binned to a pixel scale of 0'.'144 to allow for a 
comparison with the near-infrared images described below. The 
optical knots in the [SII] -image are labeled following the nota- 
tion of Ogura & Walsh ( 1992 1. A continuum subtracted image 
constructed by use of the Gunn z image shows that the optical 
knots HH136-E (IRS9 in the notation of Tamura et al. ll^97t . 
HH 136-1, and HH136-J are continuum sources with weak if 
any [SII] emission. Strong but featureless [SII] -emission oc- 
curs towards HH135. The [SII] emission at the southern bound- 
ary of the imaged region traces the bright rim of Sugitani 64. 

3.1. Near-infrared imaging of HH135/HH136 

Narrow -band images obtained in the NB212, NB164, and 
NB228 filters are shown in Figs.|5]E and|5] respectively. The 
slit locations shown in Fig. |3 are those adopted during the 
February 2004 observations and correspond to position angles 
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Fig. 4. Narrow -band image of the HH135/HH136 region ob- 
tained in the NB164 filter, which is centered on the wavelength 
of the [Fell] a'D-112 - a'^Fg/2 emission Une near 1.644 fim. The 
three double lines indicate slit locations and slit widths for the 
March 1999 spectroscopic observations, with position angles 
of PA = 34.°5, 38°, and 45°. North is up and east is left. 
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Fig. 5. Narrow-band image of the HH135/HH136 region ob- 
tained in the NB228 filter. North is up and east is left. 

of 33°, 43°, and 45°. Figure|3|includes the sHt locations for the 
March 1999 spectroscopic observations, at position angles of 
34° .5, 38°, and 45°. Observations at a position angle of 34° .5 
were repeated twice. The NB228 passband does not cover any 
strong emission lines and was used to construct the continuum- 
subtracted images in the H2 2.122 fj.m and [Fell] 1.644 fim 
lines. A scaUng factor of 0.933 that corresponds to the ratio 



of the filter widths was applied to the flux in NB228 to obtain 
the H2 image shown in Fig.0 In a similar way, the continuum- 
subtracted [Fell] image shown in Fig.|6lwas constructed using 
a scaling factor of 0.3. 

The overall morphology of the H2 emission in the 
HH135/HH136 region is that of a highly collimated bipolar 
flow with various H2 emission knots along the flow direction. 
The labeling of the infrared H2 knots of Fig.0follows the nota- 
tion of the optical knots in the [SII] image. This is possible be- 
cause in general, the overall morphology of the [SII] and the H2 
emission is very similar Exceptions occur for the optical knots 
HH136-A, C, and F, which do not have near-infrared counter- 
parts in H2 or [Fell], and which do not show significant infrared 
continuum emission either The optical knot HH136-B splits 
up into four point-like objects, B1-B4. At the NE prolonga- 
tion of knots Bl and B2, the diffuse H2 emission labeled bow- 
NE has a pronounced bow-shaped morphology that suggests a 
B2-bowNE flow direction. Knots D1-D3 (optical knot HH136- 
D) show strong H2 emission with weak underlying continua. 
Knot E is a strong continuum source (IRS-9 in the notation of 
Tamura et al. ('1997 i) without pronounced H2 emission. The 
region between infrared knots G-J2 contains a large number 
of infrared continuum sources, and very strong infrared con- 
tinuum emission occurs between knots H1-J2 (cf. Fig. |5}- In 
the NB228 continuum image shown in Fig.|5] the presence of 
a large number (30-40) of point-like continuum sources is dis- 
cerned, in agreement with the finding of Tamura et al. ('1 997> 
based on their Ks image. The slit traces shown in Figs. fTI and 
^show that most of the continuum sources that are located be- 
tween knots I1-J3 show H2 emission as well. Knots J1-J3 are 
very strong continuum sources. Because of the very strong con- 
tinuum emission between I1-J3, a proper continuum-subtracted 
H2 image for that region could not be constructed. It is noted 
that between knots J2 and d, a remarkably circular shape of 
very faint H2 emission is discerned in a false-color image, with 
knot b at the southeastern and knot a at the northeastern ex- 
tremes of this 'cavity'. 

Five infrared knots here labeled a-e are identified in 
HH135. The H2 emission of HH135 has a morphology of a 
bow with an apex near knot d and wings formed by knots b, c, 
and e. At the SW prolongation of the HH136-HH135 complex, 
the diffuse H2 emission labeled bow-SW has a pronounced 
bowshape. Faint H2 emission also occurs in the SW extreme 
of HH135. The emission has a pronounced bow morphology 
as well. The H2 continuum-subtracted image is reproduced in 
Fig. |5] The bow morphology, which cannot be demonstrated 
too convincingly in the contour plot shown here, is more ap- 
parent in a false-color H2 image. The whole extent of the col- 
limated molecular hydrogen flow is 80", which corresponds to 
about 1 parsec at a distance of 2.7 kpc. 

The morphology of the [Fell] emission is given in Fig. |6| 
The [Fell] emission is, in general, very strong, and it is associ- 
ated with molecular hydrogen emission. Exceptions occur for 
bow-NE and bow-SW, which show faint [Fell] emission only. 
The very strong [Fell] emission near knot d is displaced by a 
few arcsec towards the SW relative to its H2 counterpart. The 
displacement between the [Fell] and H2 emission provides ev- 
idence for the presence of fast, dissociative J-type shocks with 
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Fig. 6. Continuum-subtracted image showing the [Fell] 1.64 
fim emission towards HH135/HH136. 

a NE-SW flow direction, where the [Fell] emission arises in 
the fast, dissociative parts of the shocks, such as the apex of a 
bow shock, and where H2 arises in slower regions or in oblique 
shocks, such as those in the wings of a bowshock. 

The NE-extreme of the HH135/HH136 flow is shown in 
Fig. |8] together with the location of a slit at position angle 
PA=65°, which was used in February 2004 to obtain spectra 
of objects bowNE and knots Bl and B2. The image shows 
the continuum-subtracted H2 emission. Object bowNE has a 
pronounced bow shape with knots Bl and B2 located along 
a possible SW-NE flow direction. The alignment of bowNE- 
B2 differs by some 20° from the large-scale alignment of the 
HH135/HH136 complex, which suggests the presence of a sec- 
ond flow, that is not powered by IRAS 1 1 101-5829. 

The images shown in Figs. were obtained in February 
2004 at a seeing of 0'.'55. Seeing during the March 1999 imag- 
ing campaign was significantly worse and generally above 1 
arcsec. A comparison of the images obtained at both epochs 
does not show any spatial displacement of the various knots, 
which supports the assumption that HH135/HH136 is located 
at the large distance of 2.7 kpc. Thus, no effort was made 
to obtain proper motions and displacement vectors from both 
epochs. 



Fig. 7. Continuum-subtracted image showing the H2 2.122 /rni 
emission towards HH135/HH136. The strongest knots are la- 
beled. 
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Fig. 8. Morphology of the molecular hydrogen emission at the 
NE extreme of the HH135/HH136 flow. The double fines in- 
dicate tfie position of a slit and its widtfi placed at a position 
angle ofPA=65°. 
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Fig. 9. Morphology of the molecular hydrogen emission at the 
SW extreme of the HH135/HH136 flow. 



3.2. Molecular hydrogen excitation across the 
HH135/HH136 flow 

During the March 1999 and February 2004 observing cam- 
paigns, spectra of the HH135/HH136 complex were obtained 
using a total of 7 slit positions at various position angles across 
the region. The relatively large number of slit positions was 
adopted to investigate possible changes in the H2 excitation 
across the HH135/HH136 flow. The variation of the flux in the 
(1,0) S(l) line of H2 along the various slit positions is shown 
in Figs. rHO The individual emission knots, as identified in 
Fig. are marked in these slit traces. The bold lines repre- 
sent the flux variation in the continuum along the slit. It can 
be seen that many of the H2 emission line knots coincide with 
faint, underlying continuum sources. The H- and Kj-band spec- 
tra are shown in Figs.Q]and^ Note the very strong continuum 
emission that is associated with knots HI and 12. 

Tables fTI - [41 contain the unreddened line fluxes measured 
towards the various knots. From the analysis of the [Fell] emis- 
sion lines, a visual extinction of Ay = 2.7 mag is inferred 
(cf. Sect. 13. 4> and was used to obtain the population densi- 
ties A^(v'J') in the ro-vibrational levels v'J'. A(v'J'v"J") and v 
are the transition probabilities and wave numbers, respectively. 
The visual extinction obtained from the [Fell] lines is in excel- 
lent agreement with the values of Ay = 2.8 - 3.1 mag obtained 
by Ogura & Walsh ( 1992i from their optical spectroscopy. The 
population densities were used to construct excitation diagrams 
for the various knots, where values of Zn(A^(v'J')/g are plotted 
versus excitation energies ^(v'J') of the ro-vibrational levels. 
Statistical weights are g = gs(2J' + 1), where the nuclear-spin 
statistical weights gs are 1 and 3 for even and odd rotational 
levels J', respectively. The above analysis implies that the H2 
ortho/para ratio is 3, which in J-type shocks is attained at shock 



velocities above Vs ~ 15 km s or Vj > 40 km s ' for C-type 
shocks (Wilgenbus et al. l2000b . 

The excitation diagrams for the various knots are shown in 
Fig GSI Open triangles and open squares correspond to data in- 
ferred from the blue and red grisms GB and GR, respectively, 
while filled triangles and squares represent data inferred from 
the H- and K-band observations carried out at the higher spec- 
tral resolution, respectively. Straight lines fitted to the various 
population distributions show single temperature fits, with val- 
ues of 1650K for knot c, 1815K for knots b and B2, 1855K 
for knot B4, 1920K for knots Bl, B3, Dl, D3 and HI, 2030K 
for knot d and bowNE, and 2160K for knot 12. Uncertainties in 
the derived excitation temperatures are about 200 K. The tem- 
perature range among the various knots is extremely narrow, 
and except for knot c, all data are consistent with a constant 
excitation temperature of 2000 + 200K. 

In a few cases, the population densities inferred from the 
various grism settings for a given knot are not consistent with 
each other This is seen in the excitation diagrams of Fig. 1101 
where most pronounced discrepancies occur for knots Bl, Dl, 
and c. For knot Bl, population densities derived from the GB 
observations are lower by factors of a few compared to pop- 
ulation densities derived from GR and the H and Ks grisms. 
For knots Dl and c, population densities derived from the ob- 
servations with the GR grism are too low, if compared with 
the population densities derived from the GB, H, and Ks obser- 
vations. The discrepancies do not arise from a larger than as- 
sumed reddening towards these knots, which would explain the 
differences in knot Bl, but not those in Dl and c. They result 
from relatively severe slit losses during phases of poor seeing 
and from positional differences in the various grism settings, 
most pronounced for the GB and GR observations that were 
taken several years apart. However, the temperatures derived 
from the discrepant population densities agree with tempera- 
tures derived from the remaining data, for a given knot. This 
is illustrated by the dashed lines shown for knots Bl, Dl, and 
c, which are temperature fits to the discrepant population den- 
sities. Thus, discrepancies do not in any way affect the con- 
clusions derived here concerning the H2 excitation, but they 
severely hamper a more detailed analysis of the [Fell] emis- 
sion as discussed below. 

3.3. Global properties of the H2 emission in 
HH135/HH136 

As demonstrated in Fig.^| the ro-vibrational population dis- 
tribution of H2 among the various knots is remarkably similar. 
It thus makes sense to sum up the individual spectra to increase 
the signal-to-noise ratio, to search for fainter H2 emission lines 
that are not apparent in the individual spectra. This method 
simulates spectroscopic observations of HH135/HH136 with 
a very large aperture, and is thus physically meaningful. The 
sum spectra, shown in Figs. ll2tlTSl show a wealth of addi- 
tional, faint emission lines, particularly in the J-band. Some H2 
emission lines are identified explicitly in the figures. The im- 
portance of J-band observations for a proper characterization 
of the H2 excitation has been emphasized by Giannini et al. 
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Fig. 10. Excitation diagrams constructed for individual knots in HH135/HH136 with values of ln(N(v' J') / g plotted versus ex- 
citation energies E(v'J'). Open triangles and open squares correspond to data inferred from the blue and red grisms GB and 
GR, respectively, while filled triangles and squares represent the data inferred from the H- and K-band observations using the 
high-resolution grism. Straight lines correspond to single temperature fits to the population distribution, whereas dashed lines 
plotted in the diagrams of knots Bl, Dl, and c are fits to discrepant measurements (see text). 



J2002> . and full use of the additional J-band emission lines is 
made here. 

The H2 line fluxes derived from the sum spectra and 
the inferred column densities N(v'J') are given in Tables 131- 
FtI The conversion of molecular hydrogen fluxes to popula- 
tion densities was achieved via the relation /(v'J'v"J") = 
^vA(v'J'v"J")A^(v'J')- Intensities /(v'J'v"J") are inferred 
from the de-reddened fluxes F(v'J'v"J"), using the solid angles 
Q given in Sect. |2] The corresponding H2 excitation diagram 
is shown in Fig. ^2 The population distribution indicates the 
presence of a thermalized gas at = 2000 K, which contains 
a minor fraction of 0.3% of hot gas at T^.^ = 5500 K. Dots in 
the excitation diagram show the modeled population densities 
in all levels at such a gas mixture. 



Because of the wealth of [Fell] and other atomic emission 
lines in the spectra, and because of the relatively low spectral 
resolution, the confusion limit is reached in various spectral re- 
gions and care needs to be exercised when inferring line fluxes 
of individual emission lines. It thus makes sense to calculate 
a theoretical H2 spectrum and compare it with the observed 
spectrum. This allows us to judge whether the fluxes of the 
molecular hydrogen lines listed in Tables rSlTTl are severely af- 
fected from line blends. A theoretical H2 spectrum was con- 
structed for thermalized H2 in a gas mixture of 2000 K and 
5500 K with a column density ratio of 300:1, using the code of 
Gredel & Dalgarno ( 1996). The theoretical spectra, calculated 
for a visual extinction of Ay -2.1 mag, were normalized to the 
strength of the (1,0) S(l) line and convolved with a Voigt pro- 
file at the various spectral resolutions of the observed spectra. 
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Fig. 11. Excitation diagram constructed from the H2 sum spec- 
tra. Triangles, squares, and hexagrams, respectively, represent 
population distribution among rotational levels of v' = l, 2, 3, 
and 4. Dots represent a theoretical population distribution con- 
structed for a thermalized gas at a temperature of 2000K, with 
a second component added at a temperature of 5500K and a 
column density 300 times less than the first. 



The theoretical spectra are reproduced in Figs. 1121-1161 as bold 
lines. 

The comparison of the observations with the model spec- 
tra show that for single lines that appear unaffected by blends 
with atomic lines, the fit is excellent. The H2 line fluxes in the 
H- and K-band atmospheric windows are generally very well 
reproduced (cf. Figs. [Tsl and O. Note that the (1,0) S(l) line 
near 2.121 yum, as well as the a.'^Dj/2 - a''F9/2 [Fell] line near 
1.644 fj.m, show broad emission wings. In the heavily crowded 
J-band spectra, individual H2 lines that are free of blends are 
very well reproduced as well. Examples are the (2,0) S(l) and 
(3,1) S(4) pair near 1.164 //m (Fig.[l3 and die (4,2) S(5), (4,2) 
S(6), and (2,0) Q(l)-Q(4) lines at 1.21-1.25 //m (Fig.lHli. The 
very good fit of emission lines from v' =2 to v' =4 shows that the 
two-temperature fit to the H2 population distribution (Fig. II li 
is accurate, despite the flux uncertainties in many of the H2 line 
fluxes that result from line blends. 

Re-forming molecular hydrogen in the cooUng, post-shock 
regions causes a pronounced variation in the population den- 
sities of ro-vibrational levels with similar excitation tempera- 
tures, but different vibrational quantum numbers (e.g.. Flower 
et al. l2003l their Fig. 6). The effect leads to a pronounced over- 
population of rotational levels in v' =0, relative to levels in v' = 1 
of similar excitation energy. Emission lines that arise from 
v' =0 are not available in the present study. The rotational pop- 
ulation densities in v'=0 are very sensitive to C-type shocks, 
which produce very large column densities in these levels. C- 



type shocks have been used to explain the observed H2 in var- 
ious outflow sources, such as HH 99 and VLA 1623 A (Davis 
et al. 19991 or in IRAS18151-1208 (Davis et al.'2004l. Others 
have used C-type shocks with magnetic precursors (McCoey et 



al. 120041 Flower et al. l2003t . The lack of kinematical data and 
the lack of observations of H2 emission from v = makes it 
diflicult to judge whether C-type shocks are relevant for the H2 
excitation in HH135/HH136, although they must not be ruled 
out. 

For a gas mixture of warm molecular gas at 2000K plus 
a fraction of 0.3% at 5500 K, the model calculation produces 
a total H2 flux of ftot(H2) = Ev'j'vj-f (v'J'v"J") = 14.2 x 
F(1301), where f (1301) is the flux in the (1,0) S(l) line and 
where the summation is carried out over all possible emission 
fines that arise from the E2 cascade in the X^l.t electronic 
ground state of H2. The various slit positions employed here 
sample a total H2 flux of F,o, (H2) = 1.210"^^ W mT^, which 
corresponds to a total measured luminosity of L(H2) = 0.8 Lq. 
From a comparison of the H2 emission regions sampled by the 
various slit positions with the overall morphology of the H2 
emission in HH135/HH136, it is estimated that the total H2 lu- 
minosity in HH135/HH136 reaches at least 2 Lq, but possibly 
exceeds this number by factors of a few. The total H2 popu- 
lation density, for the two-component mixture assumed above, 
is Mo,(H2) = 2v-j'A?(v'J') = 77A?(1,3), or Mo,(H2) = 5.210^'^ 
cm"^. The scaling factors (14.2 and 77 in the present case) are 
dependent on the actual ro-vibrational population distribution 
and show a pronounced dependence on the H2 excitation tem- 
perature (e.g., Gredel 1994 1. 

Low-mass stars that power molecular outflows such as 
HHl 1 1 and HH212 have bolometric luminosities of a few tens 
of solar luminosities (Froebrich 2005 1, whereas the H2 lumi- 
nosities associated with such flows are of the order of 0.1 Lq. 
The accretion rates of Class and Class I objects are, pur- 
portedly, proportional to the flow luminosities (Stanke 120001 
Froebrich et al. 2003 1, and the H2 luminosities L(H2) seem to 
increase with increasing source luminosities L(bol). The very 
large luminosity of IRAS 1 1 101-5928 suggests the presence of 
a deeply embedded, intermediate- or high-mass star in forma- 
tion. Its luminosity and the luminosity of its H2 outflow define 
a data point that is consistent with the regression of L(H2) with 
L(bol) from the low-mass protostars. This result may be in- 
terpreted in terms of a significantly increased accretion rate in 
IRASl 1101-5829, that is, that the intermediate- to high-mass 
star formation witnessed by HH135/HH136 proceeds via ac- 
cretion as well, but with significantly increased accretion rates, 
relative to the accretion rates of low-mass protostars. This 
interpretation is in agreement with Brooks et al. 120031 who 
suggested that accretion in the luminous IRAS source IRAS 
16547-4247 drives the observed, well-coUimated outflow that 
extends over 1 .5 pc. Similarly, Davis et al. 120041 observed a 
collimated, parsec-scale molecular outflow from IRAS 18151- 
1208 and concluded that massive protostars drive collimated 
jets while in their earliest stages of evolution. More studies 
of intermediate- to high-mass protostars and their outflows are 
needed to understand if the outflow properties are scaled-up 
versions of the low-mass counterparts and if high-mass star for- 
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Fig. 12. Monochromatic fluxes versus wavelength of the sum spectrum including the model H2 spectrum (bold line). The sum 
spectrum is obtained from individual spectra using grism GB. A few H2 Unes are identified. 
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Fig. 13. Monochromatic fluxes versus wavelength of the sum spectrum including the model H2 spectrum (bold line). The sum 
spectrum is obtained from individual spectra using grism GB. The (4,2) S(6) transition of H2 is identified. 
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Fig. 14. Monochromatic fluxes versus wavelength of the sum spectrum including the model H2 spectrum (bold line). The sum 
spectrum is obtained from individual spectra using grism GR. The (1,0) S(l) and (1,0) S(3) Unes of H2 are identified. 
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Fig. 15. Monochromatic fluxes versus wavelength of the sum spectrum including the model H2 spectrum (bold line). The sum 
spectrum is obtained from individual spectra using the HR grism. The (1,0) S(6) and (1,0) S(7) Unes of H2 are identified. 
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Fig. 16. Monochromatic fluxes versus wavelength of the sum spectrum including the model H2 spectrum (bold line). The sum 
spectrum is obtained from individual spectra using the HR grism. The (1,0) S(0) and (1,0) S(l) lines of H2 are identified. 



mation proceeds, in general, via the same physical processes 
that govern low-mass star formation. 

3.4. The atomic emission 

The J- and H-band spectra of the various emission line knots in 
the HH135/HH136 region are dominated by a wealth of strong 
emission lines from [Fell] and from atomic species such as CI, 
Sll, and Nl. Emission from hydrogen recombination lines in the 
Brackett and Paschen series are detected as well. Optical emis- 
sion from various ions such as [Fell], [Felll], [SII], [OI], [Oil], 
[OIII], etc., and from Hel and Ha has been reported by Ogura 
& Walsh lil992j . Emission from [Fell] is generally weak in Hll 
regions and planetary nebulae, but rather strong in supernova 
remnants and in fast shocks. 

A comprehensive analysis of the [Fell] emission in Herbig- 
Haro flows has been presented by Nisini et al. (I2002t 
and by Hartigan, Raymond, & Pierson (2004 1. [Fell] emis- 
sion lines that arise from a common upper level, such as 
flie - a'*D7/2(/11.257;U m), a^Dy/a - a'*D7/2(/il.321// m), 

a^Fg/z - a'*D7/2(/11.644yu m), and a%/2 - a'*D7/2(/ll.809yu m) 
lines, can be used to determine the reddening towards the line 
emitting region. As pointed out by Hartigan et al. (2004i, the 
Einstein A-values of Quinet et al. (1996 1 are probably more 
accurate than those of Nussbaumer & Storey (1988 1 because 
the former result in a better agreement of the extinction val- 
ues determined from [Fell] lines, compared to the extinction 
determined from other methods such as the Balmer decrement. 
From the line fluxes of the sum spectra tabulated in Tables fsl 



andrgl we infer a visual extinction of Ay - 2.5 - 2.7 mag us- 
ing Quinet's oscillator strengths. This is in excellent agreement 
with the values of 2.8-3.2 mag obtained by Ogura & Walsh 
(^992 1 from their optical observations. Nussbaumer's A-values 
result in visual extinctions of 5.4-5.7 mag, which are too large 
if compared with the extinction inferred from the optical data. 
We note, however, that the good agreement between the extinc- 
tion derived here and the value of Ogura & Walsh 119921 may 
be fortuitous, as the extinction has been derived from [Fell] 
fluxes in the J- and H-band spectra, which were taken several 
years apart, and where positional differences of the various slit 
positions introduce large uncertainties in the [Fell] line fluxes 
towards some of the knots in HH135/HH136. 



Following Nisini et al. J2002> . it is in principle possible to 
estimate the electron temperature Te via the [Fell]1.257/Pa/3 
flux ratio. Our values would indicate an electron temperature 
of about Tg a: 3000 K. This value is, however, highly uncertain 
and probably wrong altogether, as pointed out by the referee. 
The Fell/PayS ratio depends on both the abundance of iron in the 
gas phase and on the ionisation fraction, two parameters that 
are not well constrained. In addition, at such low electron tem- 
peratures, the [Fell] levels of the observed lines, whose excita- 
tion energies are above 10000 K, are not populated enough to 
give rise to the observed, strong emission. The low temperature 
is also not consistent with the fast J-shocks needed to account 
for the iron presence in the gas form. The [Fell] 1 .257/[Cl]0.98 
flux ratio, compared with the [Fell] 1 .247/Pay6 flux ratio, sug- 
gests the presence of a very fast J-type shock, with shock ve- 
locities around Vj ~ 100 km s"', assuming a standard gas-phase 



12 



R. Gredel: HH135/HH136 - a luminous H2 outflow towards a high-mass protostar 



Fe abundance of Fe/H = 10 ^ The [Fell] 1.644/ 1.600 ratios are 
similar to line ratios observed in other HH objects with [Fell] 
emission. Ratios of 10-15 prevail towards most of the knots in 
HH135 and HH136 (cf. Tables[21and[3|i and indicate electron 
densities around - 4000 cm"^ or somewhat below (cf. Fig. 7 
of Nisini et al. ( I2002t '). Larger ratios of 20-30 are inferred for 
knots c and D3, and corresponding electron densities are some- 
what below rie = 3000 cm^^ Flux ratios of [Fell] 1.644/1.533 
and [Fell] 1.644/1.677 of about 10 for all the knots indicate 
electron densities of the order of = 3500 cm"^. The man- 
ifold atomic emission lines towards HH135/HH136, together 
with the low spectral resolution available here, causes a large 
number of line blends (cf. Figs. El - II6I 1. which prohibits a 
more comprehensive analysis of the atomic emission. 

4. Conclusions 

The observations presented above and the main conclusions are 
summarized as follows: 

1. The images obtained in the (1,0) S(l) line of molecular hy- 
drogen reveal the presence of a well-collimated molecular 
outflow that extents over a scale of about 1 pc. 

2. A number of the H2 emission line knots are associated with 
faint, underlying continuum emission. 

3. The ro-vibrational excitation temperatures of H2 of the 
various knots in the HH135/HH136 flow are remarkably 
constant, and are well characterized by a narrow range of 
Tex = 2000 + 200K. 

4. The molecular part of the shocked gas contains a small frac- 
tion of some 0.3% of hot H2 at a ro-vibrational excitation 
temperature of Tex = 5500 + 200K. 

5. Very strong emission lines from [Fell] occur towards var- 
ious knots in HH135/HH136, and emission from Hel, 
hydrogen recombination lines, [CI], [Sll], and [Nil] are 
present as well. The [Fell] line ratios indicate the pres- 
ence of a fast J-type shock at a speed of v, ^100 km s"'. 
Electron densities are of the order of = 3500 - 4000 
cm"^^ and electron temperatures are ^ 3000 K. A more 
comprehensive analysis of the atomic emission is hampered 
by the low spectral resolution and the occurrence of a very 
large number of line blends. 

6. The ro-vibrational population distribution of H2, together 
with the presence of strong [Fell] emission, which is spa- 
tially displaced from the H2 emission, indicate that the 
emission lines arise form in the cooling regions of fast, dis- 
sociative J-type shocks, where the [Fell] emission traces 
the fast, dissociative parts of the shocks and where the H2 
emission emerges from regions of oblique shocks where the 
shock speeds are lower The presence of emission from re- 
forming H2 molecules in the cooling post-shock region is 
not ruled out. 

7. The large H2 luminosity of 2 Lq suggests that the 
intermediate- to high-mass protostar that powers the 
HH135/HH136 outflow forms via a significantly increased 
accretion rate, compared to the accretion rates of low-mass 
protostars. 
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Fig. .1. Spatial variation of the (1,0) S(l) H2 emission line flux 
along the sUt for the three positions PA=34.°5, 38°, and 45° as 
inferred from the March 1999 observations. Bold lines show 
the variation of the continuum emission. 



Fig. .3. H-band spectra of individual knots in the HH135 flow. 
Spectra are shifted by units of 1000, 2000, 3000, 4000, 5000, 
6000, and 8000 for knots HH135-B1, B4, Dl, D3, HI, 12, and 
HH136-d, respectively, along the ordinate. 
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Fig. .4. K-band spectra of individual knots in the HH135 flow. 

The ordinate is valid for bowNE only. Spectra are shifted by 
units of 1000, 2000, 3000, 4000, 5000, 6000, and 8000 for 
knots HH135-B1, B4, Dl, D3, HI, 12, and HH136-d, respec- 
tively, along the ordinate. 



Fig. .2. Spatial variation of the (1,0) S(l) H2 emission line flux 
along flie sHt at the three positions PA=33°, 45°, and 43° as 
inferred from the February 2004 observations. Bold lines show 
the variation of the continuum emission. 
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Table .5. H2 line fluxes and column densities in ro- vibrational 
levels v'J', as measured from the low-resolution sum spectinm 
obtained with the blue grism GB. 



Line transition 


Wavelength 


Flux 


A'(v'J') 






10"^" W m"- 


lO^'cm-- 


(2,0)S( 9) 


1.054 


156 


7(1) 


(2,0)S( 7) 


1.064 


128 


5(1) 


(2,0)S( 6) 


1.073 


141 


6(1) 


(2,0)S(4) 


1.100 


82 


4(1) 


(2,0)S( 3) 


1.117 


379 


19(10) 


(2,0)S( 2) 


1.138 


225 


14(7) 


(3,1)S(5) 


1.152 


188 


4(2) 


(2,0)S( 1) 


1.162 


288 


22(11) 


(3,1)S(4) 


1.167 


55 


1(0.5) 


(3,1)8(3)° 


1.186 


170 


4(1) 


(4,2)S( 9) 


1.196 


62 


1(0.5) 


(4,2)S( 8) 


1.199 


59 


1(0.5) 


(4,2)8(7) 


1.205 


67 


1(0.5) 


(4,2)8( 2) 


1.285 


56 


1(0.5) 


(4,2)8( 6) 


1.214 


27 


0.4(0.2) 


(4,2)8( 5) 


1.226 


69 


1(0.5) 


(3,1)8( 1) 


1.233 


156 


5(1) 


(2,0)Q( 2) 


1.242 


95 


11(2) 


(2,0)Q( 3) 


1.247 


465 


57(11) 


(2,0)Q( 4) 


1.255 


251 


32(6) 


(2,0)Q( 5) 


1.264 


420 


55(11) 


(2,0)Q(7) 


1.287 


114 


16(3) 


(3,1)Q( 1) 


1.314 


138 


5(1) 


(4,2)8( 1) 


1.312 


125 


3(1) 



Possible blend with [PI]. 



Table .6. H2 line fluxes and column densities in ro- vibrational 
levels v'J' , as measured from the low-resolution sum spectrum 
obtained with the red grism GR. 



Line transition 


Wavelength 


Flux 


iV(v'J') 




fim 


10-2"W m-^ 


10''*cm-2 


(1,0)8(9) 


1.688 


510 


65(13) 


(1,0)S( 8) 


1.715 


700 


65(13) 


(1,0)8(7) 


1.748 


2798 


207(41) 


(1,0)S( 3) 


1.958 


17948 


1056(528) 


(1,0)8(2) 


2.034 


3016 


195(97) 


(2,1)8(3) 


2.074 


897 


41(20) 


(1,0)8( 1) 


2.122 


8657 


669(134) 


(2,1)8(2) 


2.154 


416 


20(4) 


(1,0)8(0) 


2.223 


2599 


289(58) 


(2,1)S( 1) 


2.248 


521 


30(6) 


(1,0)Q( 1) 


2.407 


10675 


757(454) 


(1,0)Q( 2) 


2.413 


6122 


616(370) 


(1,0)Q( 3) 


2.424 


8936 


984(590) 


(1,0)Q(4) 


2.437 


4419 


514(308) 
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Table .1. Individual measurements of emission lines towards the various knots in HH135 and HH136, using the low resolution 
spectra obtained with grism GB. 



Wavelength (pm) Line flux(10-20W m"^) 

PA33 PA43 PA65 





D 1 


DA 




DA 


Ul 


ML 


e 


a 


D 1 


DO 
D/ 


0.953 


0.4 


1.4 


1.0 


0.3 


0.6 


— 


— 


3.5 


0.6 


0.6 


0.983 


0.2 


0.3 


0.3 


0.1 


0.3 


0.1 


0.3 


0.7 


0.3 


0.3 


0.986 


0.5 


0.9 


0.8 


0.4 


0.6 


0.2 


0.8 


1.9 


0.7 


0.6 


1.006 


0.2 


0.5 


0.6 


0.2 


0.4 


0.2 


0.1 


1.3 


0.4 


0.3 


1.0293 


0.5 


1.7 


1.8 


0.4 


1.8 


0.1 


0.3 


4.6 


0.7 


0.6 


1.033 


1.3 


3.8 


4.2 


1.1 


2.6 


0.4 


0.8 


10.9 


1.6 


1.5 


1.038 


0.3 


0.8 


0.8 


0.2 


0.8 


0.1 


0.1 


2.0 


0.3 


0.4 


1.041 


0.5 


1.2 


2.1 


0.2 


1.2 


0.2 


0.2 


4.9 


0.4 


0.4 


1.053 


0.1 


_ 


0.4 


0.1 


_ 


_ 


0.1 


0.3 


_ 


_ 


1.064 


0.2 


0.3 


0.2 


— 


0.1 


— 


0.2 


0.5 


- 


0.2 


1.072 


0.2 


0.4 


0.3 


- 


- 


- 


0.2 


0.9 


0.3 


0.2 


1.084 


3.3 


10.2 


8.7 


3.4 


6.8 


1.1 


1.8 


22.7 


4.4 


3.9 


1.095 


0.3 


1.0 


1.1 


0.5 


1.0 


0.4 


0.5 


2.4 


0.4 


0.5 


1.163 


0.2 


0.4 


0.4 


0.1 


0.3 


0.1 




0.6 


0.3 


0.3 


1.189 


0.4 


1.1 


1.2 


0.4 


0.8 


0.2 


0.3 


2.5 


0.8 


0.6 


1.233 


0.1 


0.2 


0.2 


0.1 


0.2 


0.1 


0.3 


0.5 


0.2 


0.1 


1.239 


0.1 


0.3 


0.3 


0.2 


0.2 


0.1 


0.4 


0.4 


0.3 


0.3 


1.248 


0.3 


0.7 


0.9 


0.3 


0.9 


0.2 


0.3 


1.7 


0.6 


0.5 


1.257 


8.5 


20.0 


22.4 


7.7 


17.1 


3.8 


5.9 


39.5 


11.1 


9.8 


1.271 


0.4 


1.2 


1.5 


0.4 


1.1 


0.1 


0.1 


2.6 


0.5 


0.4 


1.279 


0.7 


2.0 


2.4 


0.9 


2.2 


0.3 


0.2 


4.4 


0.9 


0.7 


1.282 


0.7 


2.4 


3.8 


0.6 


2.8 


1.4 


1.0 


6.3 


1.2 


1.2 


1.295 


1.1 


2.9 


3.7 


1.0 


2.8 


0.5 


0.5 


6.6 


1.4 


1.1 


1.298 


0.3 


0.8 


0.9 


0.2 


0.7 


0.1 




1.7 


0.3 


0.3 


1.321 


2.5 


5.7 


7.3 


2.2 


5.6 


1.3 


1.8 


12.0 


3.4 


2.9 


1.328 


0.7 


1.8 


2.1 


0.5 


1.6 


0.2 


0.3 


3.7 


0.8 


0.5 


1.534 


1.7 


4.2 


6.6 


1.4 


4.9 


1.1 


0.8 


8.4 


2.0 


1.6 


1.600 


1.2 


3.1 


5.5 


1.0 


4.0 


0.8 


0.6 




1.5 


1.3 



Flux uncertainties are 10% in general, but increase to about 50% for I.IO <A< 1.18jum. 
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Tkble .2. Individual K-band measurements of emission lines towards the various knots in HH135 and HH136, using the medium 
resolution spectra obtained with grism HR. 



Wavelength (jum) Line flux(10-20W m"^) 

PA65 PA45 PA43 PA33 





B2 


fil 


bowNE 


d 


c 


01 


D3 


S3 


bowNE 


h 


HI 


12 


01 


S4 


SI 


2.034 


14.2 


8.4 


2.9 


17.0 


8.4 


21.3 


6.6 


6.5 


2.9 


7.6 


17.9 




37.1 


12.5 


10.7 


2.046 


0.9 






1.0 


0.2 


3.6 




0.3 










7.8 


1.4 




2.074 


2.6 


2.1 


0.6 


2.9 


0.8 


5.2 


1.4 


1.1 


0.6 


1.6 


2.9 




7.3 


2.9 


2.3 


2.122 


39.5 


26.3 


6.8 


47.9 


21.0 


59.2 


19.0 


19.5 


7.3 


22.7 


48.4 


19.6 


103.3 


36.3 


29.9 


2.154 


1.4 


1.0 


0.4 


1.7 


0.8 


2.5 


0.8 


0.5 


0.3 


1.7 


1.0 




4.0 


1.5 


1.1 


2.166 


0.6 


0.4 




0.2 


0.5 


3.6 


1.4 


0.7 




1.2 




1.5 


7.1 


1.6 




2.201 


0.7 


0.5 


0.2 


0.5 


0.4 


1.2 


0.3 


0.3 


0.4 




0.4 




1.5 


0.9 


0.6 


2.223 


11.1 


7.4 


1.7 


12.5 


5.7 


20.8 


6.6 


6.5 


1.6 


6.6 


13.8 


5.7 


37.7 


12.4 


8.0 


2.247 


3.4 


2.4 


0.6 


4.5 


3.1 


6.6 


2.4 


1.9 


1.0 


2.1 


4.2 


2.3 


10.5 


3.5 


2.9 



Flux uncertainties are 10% in general, but increase to about 50% for A < 2.05/<m. 



Table .3. Individual H-band measurements of emission Unes towards the various knots in HH135 and HH136, using the medium 
resolution spectra obtained with grism HR. 



Wavelength (jum) Line flux(10-^"W m"^) 







PA33 






PA43 






PA33 




PA45 


PA65 


S1S4 


DID3 


HI 


S3 


Dl 


HI 


e 


S4 


01 


HI 


c 


d 


SI S2 


1.534 


12.5 


18.7 


1.9 


3.6 


9.8 


3.4 


1.9 


3.0 


8.3 


3.1 


2.8 


23.5 


2.0 6.2 


1.600 


9.3 


12.7 


1.1 


2.3 


7.4 


1.9 


1.2 


1.4 


5.3 


2.1 


1.5 


15.7 


1.7 4.8 


1.644 


99.3 


159.4 


24.4 


32.2 


71.0 


30.7 


20.9 


29.2 


64.9 


29.5 


26.7 


155.4 


23.4 65.3 


1.664 


6.4 


7.8 


0.9 


1.9 


5.3 


1.4 


0.7 


0.6 


3.1 


1.3 


1.3 


10.9 


1.2 4.3 


1.677 


10.5 


17.3 


2.0 


3.2 


8.9 


2.8 


1.5 


2.0 


7.3 


2.6 


2.0 


19.7 


2.1 6.8 


1.688 


2.0 


1.5 


1.3 




2.0 


1.5 


1.1 


0.5 


0.6 


1.3 


1.6 


3.5 


0.3 1.7 


1.712 


2.8 


3.5 


0.7 




2.2 


1.0 


0.1 


0.2 


1.6 


0.6 


0.6 


6.1 


0.4 1.3 


1.715 


1.2 


0.6 






0.6 


0.6 


0.7 




0.3 


0.5 


1.0 


2.2 


- 0.7 


1.736 
























2.4 




1.745 


2.3 


4.5 






2.5 


0.6 


0.3 




2.4 


0.8 




6.9 


0.6 2.1 


1.748 


11.0 


12.6 


3.6 


2.2 


5.4 


5.5 


5.6 


1.3 


4.3 


5.5 


6.5 


14.2 


2.0 6.9 


1.788 


5.9 


6.8 






2.8 


3.6 


3.8 




2.3 


6.0 


4.1 


7.2 


1.2 4.2 


1.798 


4.8 


5.5 






3.8 








3.2 


3.3 


0.5 


5.7 


- 3.4 


1.800 


1.8 


0.6 






3.8 








0.8 






4.1 


- 2.6 


1.810 


5.6 


9.5 






5.6 








4.2 






8.1 


- 4.3 



Flux uncertainties are 10% in general. 
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Table .4. Individual H-band measurements of emission lines towards the various knots in HH135 and HH136, using the low 
resolution spectra obtained with grism GR. 



Wavelength (jxm) Line flux(10-20W m"^) 

PA38 PA45 PA65 PA34.5 PA34.5 





Dl 


Bl 


BA 


01 


d 


c 


d 


howNE 


B2 


D3 


m 


b 


03 


m 


b 


1 


H.J 


9 s 


8 n 

O.U 


1 1 A 


90 9 
zu.z 


U.J 


1 s n 

1 J.U 


n 9 
u.z 


4 s 
4. J 


n Q 
yt.y 






1 9 
1 .z 


J.U 


n 1 

U. 1 


1 /^nn 
1 .ouu 


J.U 


1 S 
i .0 


J.J 


8 7 
0. / 


14 9 


U.J 


1 1 1 




J.O 


U.J 






n 4 
u.4 


1 Q 
i .y 




1 f\dA 
^ .Oft 


77 ^ 


1 8 (\ 


54 n 


Q1 7 
y L, 1 


119 9 

1 IZ.Z 


8 4 


89 8 
oz.o 


9 4 


JU. / 


7. J 


fl 8 
u.o 


J.J 


19 1 
iz. 1 


9*1 n 

Zj.U 


9 9 
z.z 




1.0 


n 8 

u.o 


J.J 


J.U 


7 4 




1 

0. 1 




1 7 
1. / 


U.J 






n 4 

u.4 


n 

U.7 




i.O /o 


J.J 


1 a 
1.0 


1 

0. 1 


1 n 7 

lU. / 


1 S 9 

1 J.Z 


u.o 


1 1 4 


n 9 
u.z 


^ 8 
J.O 








n 7 

U. / 


9 "X 

Z.J 




1719 


U.o 




1 4 


9 ^ 

Z.J 


4 ^ 
^. J 


n 1 
u. 1 


J.VJ 


n 9 
u.z 


1 n 
1 .u 


n 9 
u.z 






n 9 
u.z 


n 7 
u. / 




1 7zLS 


^ 7 


J.J 


7 Q 
/ .y 


1 9 1 


1 8 ^ 
10. J 


9 n 


14 Q 

14.7 


1 9 
1 .z 


4 7 


1 Q 
1 .y 






1 f\ 

i.O 


A 4 
0.4 


1 9 
i .z 


1 78Q 
1. 1 oy 


1 A 


1 7 
1. / 


^ 9 
J.Z 


S 8 

J.O 


J.J 


1 4 


D.J 


n 7 
u. / 


1 8 
1.0 


n 8 

u.o 






9 9 
z.z 


1 

j.i 


1 Q 


1 708 




1 « 

i.j 


4 8 


1 9 

1 J.Z 


1 9 

1 J.Z 


1 4 


4 


n 4 

U.4 




J. 7 


9 ^ 

Z.J 






J.O 


9 

J.Z 




1 Sin 

i . 1 u 




14 Q 


^4 a 


76 Q 


8S 7 

OJ. / 


u.u 


79 ^ 

/Z.J 


9 9 
z.z 


9fi 7 
zu. / 


1 n 8 

lU.O 






1 S S 

1 J. J 


17 S 
1 / . J 




1 


99 n 


9/1 ^ 

24. J 


^7 Q 


7/1 A 
/4.0 


0/1 ^ 


A 
0.4 


1 o 


7 /I 
/.4 


99 ^ 
ZZ.J 


7 7 
I.I 






0.0 


jO.O 




2.034 


4.7 


5.5 


10.3 


19.8 


18.4 


2.8 


17.2 


1.7 


7.2 


2.9 


5.7 


2.9 


6.4 


14.2 


5.8 


2.047 


0.7 












5.6 




0.4 








0.7 






2.074 


0.7 


0.9 


2.0 


2.8 


2.6 




3.4 


0.3 


1.3 


0.2 








5.5 




2.122 


12.1 


14.5 


29.9 


53.0 


48.5 


7.2 


47.1 


3.8 


19.1 


7.1 


10.4 


7.3 


12.4 


27.7 


12.7 


2.155 


0.5 


0.6 


1.2 


1.9 


3.2 


0.7 


1.0 


0.3 


0.6 


0.2 






0.5 


1.7 




2.167 


1.1 




0.8 


3.3 


4.1 


0.3 


1.4 






0.4 






0.7 


2.9 




2.202 


0.3 


0.3 


0.7 


1.1 


1.3 


0.3 


0.8 


0.2 
















2.224 


4.3 


4.5 


9.7 


18.4 


17.6 


2.1 


14.9 


1.2 


5.8 


2.4 


4.0 


3.2 


4.9 


8.5 


3.7 


2.248 


1.3 


1.6 


3.0 


5.6 


6.0 


1.0 


4.9 


0.6 


1.5 


1.0 


1.9 


1.6 


2.0 


3.7 


1.3 


2.355 


1.3 


0.5 


0.7 


1.9 


2.2 


0.4 


1.7 


0.5 


0.5 


0.6 












2.406 


13.5 


15.4 


31.1 


58.9 


50.9 


7.7 


53.3 


4.1 


21.6 


8.9 


13.1 


7.9 


16.4 


26.7 


15.3 


2.413 


5.8 


8.7 


15.2 


36.5 


25.8 


4.0 


26.8 


2.0 


10.0 


4.1 


6.3 


4.3 


7.9 


12.3 


7.6 


2.423 


14.0 


15.4 


31.6 


59.2 


50.9 


7.2 


53.1 


4.1 


21.0 


7.2 


12.4 


9.1 


13.8 


29.2 


16.7 


2.437 




5.0 


9.0 


16.9 


15.3 


3.3 


17.6 


1.5 


6.2 


2.4 












2.453 






1.3 








4.5 




2.3 














2.475 




2.1 


3.8 


12.9 


4.9 




4.9 


















2.498 




3.7 


9.7 


25.4 


12.5 




12.5 



















Flux uncertainties are 10% in general, but increase to about 50% for 1.9 < X< 2.05jum and for X > 2.43//m. 
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List of Objects 

HH135 [ on page 2 
HHOef on page 2 
IRASl 1101-5921 on page 2 



Table .7. H2 line fluxes measured from the medium-resolution 
H- and K-band sum spectra obtained with the HR grism in or- 
ders 4 and 3, respectively. 



Line transition 


Wavelength 


Flux 


iV(v'J') 




yum 


10 W m 


inl4 -2 

10 cm 


(1,0)S( 9) 


1.688 


652 


83(17) 


(1,0)S( 8) 


1.715 


278 


26(5) 


(1,0)S(7) 


1.748 


2328 


172(34) 


(1,0)S(6) 


1.788 


1504 


96(19) 


(1,0)S( 2) 


2.034 


2977 


192(96) 


(3,2)S( 5) 


2.066 


198 


11(6) 


(2,1)S( 3) 


2.074 


623 


28(14) 


(1,0)S( 1) 


2.122 


8763 


677(135) 


(2,1)S(2) 


2.154 


371 


18(4) 


(3,2)S( 3) 


2.201 


159 


8(2) 


(1,0)S(0) 


2.223 


2901 


323(65) 


(2,1)S( 1) 


2.248 


970 


55(11) 
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Table .8. Line fluxes measured for atomic transitions using the 
low-resolution spectrum obtained with grism GB. 



Transition 


Wavelength 


Flux 




Jim 








QICl 

ooy 


reli: 


V.yjj 


1 T? 

ill 




V.yoj 


1 'TO 

1 /V 


L>2 - r2 


u.yoj 


0/0 


XJ T "2 

xl /—J 


I.UUj 


j4U 


7 


1 OIQ 


uz. 


2d 2 CTT 


1 noo 
L.vZy 


luiy 


2P 2 Pi CTT 


1 


9AQ7 
Z4V / 


2d 2 r-\ CTT 
1^1/2 - L)5/2 


i.KJJO 


All 
4ij 


2d 2 7~\ MT 
J^l/2 - ^5/2 JNl 




1 m 
luiy 


riei 


i.Voj 




rl 0—:) 




OlO 


5d 5 c r\i 




lyj 


7 


1 1 "2 /I 
1.1 j4 


lyz 


7 


1 AAA 
1. l'f'4- 


Dy 


DTT 


1 1 /IT 
1.14/ 


1 1 A 
1 lo 


7 


1 1 ztQ 

1.147 


:?y 


7 


1.1 /O 


JO 


a U9/2 - a U5/2 


1101 

i.iyi 


14j 


TTaTT 

reli 






a JJ9/2 - a JJ7/2 


l.Zj / 


111 'XA 
1 1 1 j4 


CaTT 






«■ ^^1/2 a -I--'i/2 


1.271 


684 


Fell 






a^D3/2 - a'^D3/2 


1.279 


1132 


Fell 






H5-3 


1.282 


1645 


a^D5/2 - a4D5/2 


1.294 


1650 


Fell 






a^D7/2 - a'^D7/2 


1.321 


3595 


Fell 






a^D3/2 - a4D5/2 


1.328 


1083 


Fen 






a'^F9/2 - a'*D5/2 


1.533 


2439 


Fen 






a^F7/2 - a'^D3/2 


1.600 


1866 



Fen 



Table .9. Wavelengths and hne fluxes of atomic transitions, de- 
rived from spectra obtained with grisms GR and HR in orders 
4 and 3, respectively 



Transition 


Wavelength 


flux 

10-2"W m-2 


a^Fg/z - a^Ds/z 


1.534 


2459 


Fell 






a F7/2 - a D3/2 


1.600 


1380 


Fell 






H 12^ 


1.641 


200 


a*F9/2 - a*D7/2 


1.644 


13521 


Fell 






a'*F5/2 - a'*Di/2 


1.664 


563 


Fell 






a'*F7/2 - a'*D5/2 


1.677 


1512 


Fell 






H 1 1-4 


1.681 


568 


'Po D Hel 


1.700 


176 


T T 1 A 


1.737 


972 


H 7-4 


2.166 


1191 


erF9/2 - a'*D5/2 


1.533 


2637 


Fell 






a F7/2 - a*D3/2 


1.600 


2130 


Fell 






a*F9/2 - a''D7/2 


1.644 


13830 


Fell 






a''F5/2 - a''Di/2 


1.664 


1356 


Fell 






a''F5/2 - a''D5/2 


1.677 


2597 


Fell 






H 10-4 


1.736 


620 


a''F3/2 - a''Di/2 


1.745 


1295 


Fell 






a''F3/2 - a''D3/2 


1.797 


1017 


Fell 






a''F5/2 - a''D5/2 


1.800 


1301 


Fell 






a''F7/2 - a''D7/2 


1.809 


2173 


Fell 






H7-4 


2.166 


613 



